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\  1.  INTRODUCTION 

In  an  Influential  and  controversial  paper.  Raup  and  Sepkoski  £2^  defined  an 
event  of  mass  extinction  to  have  occurred  in  any  time  period  (of  roughly  6.4 
million  years)  for  which  the  data  value  for  that  time  period  (equal  to  the 
proportion  of  the  families  existing  at  the  beginning  of  that  period  that  went 
extinct  during  the  period)  exceeded  that  of  its  inmedlate  neighbors.  In' this 
note  wj  analyze^ the  occurrence  of  such  events  when  the  data  are  randomly 
generated  from  a  continuous  distribution. 
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Let  Xj.Xg. —  be  a  sequence  of  Independent  and  identically  distributed 
continuous  random  variables,  and  define  an  event  to  occur  at  time  period  j  if 


Xj-1  <  Xj  >  Xj+1  ,j  *  2-  Let 


j  1  if  an  event  occurs  at  time  j 

j  0  otherwise 


and  set 


N(n)  =21 
J=2  J 


In  words,  N(n)  is  the  number  of  events  by  time  n. 


L:  With  probability  1. 


lim  N(n)/n  =  1/3 


Proof :  This  follows  from  the  ergodic  theorem  since  {Ij,  j  ^  2}  is  a 


stationary  ergodic  sequence  and  E[Ij]  *  1/3.  A  simpler  proof  is  obtained  by 


noting  that  { I^ • ,Ig»...},  {I^.Ig.Ig,...},  { 1^ , 1^ ■ * • • }  are  each  lid 


sequences  whose  average  value  converges,  by  the  strong  law  of  large  numbers. 


to  1/3.  QED 


The  swan  and  variance  of  N(n)  are  easily  obtained: 


Accesion  For 


E[N(n)]  =  2  E[I.]  =  (n-l)/3 
J=2  J 
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HA 


Var[N(n)]  =  2  Var(I .)  +  2  2  2  Cov(I . . I .) 

J=2  3  i<j  J 


Var(I j)  = 

Cov(Ii * Ii+1)  = 

Cov(Ii.I1+2)  = 


-  ECI^  E[Ii+1]  =  -1/9 

P{Xi_1  <  Xl  >  Xi+1  <  X1+2  >  X1+3}  -  1/9 

S  S  S  S  S  dx.  dx,,  dx_  dx .  dx_  -1/9  =  1/45 

X1<VX3<X4>V 

o<Xl<i 


Cov(Il«Ii+j)  =0.  i  13 


we  see  from  the  above  that,  for  n£3. 


Var[N(n)]  *  2(n-l)/9  -2(n-2)/9  +2(n-3)/45 
=  (2n  +  4)/45 


It  follows  from  the  above  and  the  central  limit  theorem  for  stationary  ergodic 
processes  that 


N(n)  -  (n-1 )/3 


N(0.1) 


'/T2n+4)?45" 


t.t! 


V  j" 


W1 


V 


Terminology  If  X^_j  <  >  X^+^  call  X^+^  the  end-of-peak  value. 

3.  INTEREVENT  TIMES 

Let  Tj.  i£l,  denote  the  time  between  the  i**1  and  (i+l)st  event  -  for 

instance,  if  I2=I^=1  then  T^=2.  Even  though  the  Tj  are  neither  independent  nor 

identically  distributed  we  will  show  that  the  proportion  of  them  that  equal  j 

will,  with  probability  1,  converge  to  a  constant  value  which  we  shall  call  p  . 

J 

To  determine  the  values  p^,  call  an  event  a  J-event  if  the  event  preceding  it 

occurred  exactly  J  time  periods  earlier,  and  note  that  by  the  same  argument  as 

used  in  Proposition  1.  the  rate  at  which  J-events  occur,  call  it  r(j),  will 

equal  the  limiting  probability  of  a  j-event  at  time  pteriod  n.  Let 

R(k)  =  2  r( j)  . 

J>k 

That  is.  R(k)  denotes  the  rate  at  which  interevent  times  of  length  greater 

than  k  occur.  Now.  an  lnterevent  time  will  exceed  k  if.  for  some  r=0,l . k  , 

the  next  r  values  after  an  end  of  peak  are  each  less,  and  the  next  k-r  values 
are  each  greater,  than  their  preceding  values.  Therefore,  for  k£2. 


R(k)  =  P{X_2<X_1>XQ<X1<...<Xk}  +  P{X_2<X_1>X0>X1>...>Xk} 
k-1 

+  2  p<x_2<x_1>x0>x1>..>xrar+1<...<xk} 

r=l 


=  f  n-*>k  n-*2>  hv 

Ok!  2 


+ 


1  k+2  k-1  1  x 

S  - _ dx  +  2  J*  J 

0  (k+1) !  r=l  0  0 


l^Zirx 

r ! 


/ i  xk-r 
^  ^  dy  dx 

(k-r)! 
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where  the  last  equality  obtains  from  conditioning  the  first  probability  on  the 
value  of  XQ.  the  second  on  X_j.  and  the  remaining  on  both  X_^  and  Xf 


As  the  rate  at  which  events  occur  is  1/3  it  follows  that  the  proportion 
of  events  that  are  j-events  is  3r(J).  That  is, 

Pj  =  3  r(j). 

A  direct  Integration  of  the  above  yields,  upon  using  the  identity 

k 

2  r(j)  =  1  -  R(k)  . 

J=2 

that 

P2  =  2/5 

P3=  1/3 

P4  =  6/35 
p5  =  1/15 
P6  =  .02116401 

and  so  on.  The  quantities  Pj  also  represent  the  limiting  distribution  of  the 
Tj.  A  simple,  though  tedious,  summing  of  the  values  R(k)  yields  the  (already 
known)  result  that  the  mean  of  this  limiting  distribution  is  3. 

Remark:  The  above  could  also  have  been  obtained  by  consideration  of  the  Markov 

Chain  of  successive  end-of-peak  values.  The  limiting  value  of  Tj  could  then  be 

computed  by  conditioning  on  X,  the  limiting  end-of-peak  value  (whose  density 

o 

is  given  by  f(x)=3(l-x  )/2,  0<x<l),  and  then  computing  the  probability 
density  of  the  time  between  peaks  when  the  initial  one  has  an  end-of-peak 
value  X. 
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